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SUMMARY 


A derivation of the form of the incompressible turbulent skin- 
friction law for an insulated flat plate is made in such a way that it 
may be extended to compressible flows. The ratio of compressible to 
incompressible skin friction is obtained, and the results are shown to 
be in agreement with existing experimental results. 


INTRODUCTION 


The magnitude of the skin- friction drag encountered by a flat plate 
immersed in a fluid at Reynolds numbers large enough to insure turbulent 
flow has been one of the basic problems of aerodynamics. From the theo- 
retical approaches of Prandtl and von Karman (for a resume, see refer- 
ence 1) and the experimental work of numerous investigators, principally 
Nikuradse and Ludwieg and Tillmann (see references 2 and 3) , much has 
been learned of the forms of the turbulent boundary- layer velocity pro- 
file and the skin friction associated with these forms in incompressible 
flows, although the exact mechanisms involved are still not completely 
understood. 

Recently, the magnitude of the turbulent skin friction on a flat 
plate at high Mach numbers has become of great interest, and several 
papers have been written on this subject presenting both theoretical 
treatments of the problem and the results of skin-friction measurements 
at Mach numbers between 1.5 and 3*0. (See references ^ to 8.) 

The agreement between these theories and the experimental data that 
exist is, in general, satisfactory. The status of the problem, however, 
is such that a simple physical approach to the extension of the incom- 
pressible skin- friction laws to the compressible case would seem desir- 
able. The purpose of this paper is to present such a simple physical 
picture. 
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SYMBOIS 



roughness height — 

constant relating mixing length to y 
mixing length 
constants 
Mach number 

definitive ratio of total shear stress to laminar shear stress 


Reynolds number 


Reynolds number 


Reynolds number 



absolute temperature 
velocity in x-direction 
velocity in y-direction 
distance -along surface 
distance normal to surface 
ratio of-specific heats 


8 


boundary-layer thickness 
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p viscosity 

V kinematic viscosity 

p density 

t total shear stress 

Subscripts: 

L conditions at edge of laminar sublayer 

o free-stream conditions 

lam stress produced by laminar action alone 

turb stress produced by turbulent action alone 

w wall conditions 

THE INCOMPRESSIBLE SKIN-FRICTION LAW 


If the shearing stress due to turbulence 
velocity gradient is assumed to be given by 


in a region subject to a 


T turb 


P l 


2 du 

dy 


du 

dy 


which is the mixing-length formula of Prandtl (see reference 1, p. 130), 
the equation of motion to be satisfied by the boundary layer on a flat 
plate is 


u 
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<3u _ l Bt _ 1 5 
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Many investigators have established that for most purposes 

\ 



( 2 ) 
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is an adequate representation of the "boundary- layer velocity profile 
outside a very thin laminar sublayer. For a flat plate in the absence 
of pressure gradient 'the value of n is approximately 7 . It is also 
well-known that the value of the mixing length in the boundary layer 
near the laminar sublayer is given by 


l = ky 


(3) 


where k is approximately 0.4. 

When the preceding facts are used, the ratio of the total stress 
to the laminar stress at points in the boundary layer near the wall 
may be found in the following manner: 

Since 


T 

T lam 



the insertion of 1 from equation ( 3 ) and du/dy obtained from equa- 
tion ( 2 ), namely, 


1 -n 

Su = u oF n 

dy 1 

n5 n 


yields 


T 

T lam 


n +1 

u 0 k 2 y h 

1 + — 

nvS n 


(4) 


This equation is consistent with the idea that the total stress in 
the boundary layer becomes principally laminar as the surface of the 
plate is approached. The extent of this laminar sublayer may be reckoned 
by computing the value of y where the ratio of the total stress to the 
laminar stress is given by some definitive ratio. When this definitive 
ratio of total stress to laminar stress is assumed to be r, then from 
equation (4) 
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u Q k28 L * 
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nvS n 


and the relative thickness of the laminar sublayer is 


£l _ fn 

6 


n(r - l) v _ | n+1 


The velocity in the boundary layer at this point may be computed from 
equation (2) and is found to be 


_ n(r - 1) v~] : 

u o _ k 2 Uo5| 


In the region 0 < y % Sj, the flow is principally laminar and. 


since -i- = p. 


= 0 at the wall as may be seen from equation (l). 


the velocity profile is essentially a straight line in this region. 

The boundary layer is now assumed, for the purpose of determining the 
wall stress, to be sharply divided into a turbulent region having a 
power-law velocity profile and a laminar region having a straight-line 
velocity profile as shown in figure 1 by the solid lines. The real 
state of affairs is indicated approximately by the dashed line. The 
stress in the laminar layer just below the assumed intersection, the 

point where 6 = 5 l and u = u-^, is now given by which will 

St 

be taken as pu-^^6^. Since = 0 at the wall or the stress in the 
laminar sublayer is approximately constant. 


_ M- u L _ ^ u o [~n 
S L " 6 |_ 


(r - 1) 
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The skin-friction coefficient is 



For n = 7> equation (9) becomes 


(9) 


Cf- = 


2 


~ k2 
T(r - 1) 



1 


= Constant R 



(10) 


Equation (10) is one form of the skin-friction law that has been 
shown to be very close to the true state of~ affairs for the turbulent 
boundary layer. The value of-the constant found experimentally is 0.045. 

(See reference 1, p. 147. ) mi, “ — — — 


The value of the ratio 
yield this result ^-is 0.0444. 


l' 


which will 


The assumption of a definitive ratio of stresses is tantamount to 
assuming that the Reynolds number u^ which is characteristic of 

the laminar sublayer, is a constant. This Reynolds number formed from 
the values of 6^ and u^ from equations (6) and (7) is 


„ u L 5 L n(r - 1) 
% = — = — J2— 


Putting n = 7 


and 


k2 

r - 1 


0.0444 into equation (ll) gives 


( 11 ) 


R L = 158 


It is interesting to consider the significance of this number. It 
may well be close to the critical Reynolds number of the straight-line 
velocity profile in the laminar sublayer. Below Such a Reynolds number 
all disturbances would be damped; whereas above such a Reynolds number 
disturbances would feed on energy from the steady motion. This condition 
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would define the laminar sublayer thickness and also a mechanism by 
which the turbulent boundary layer might receive energy in order to 
continue its existence. 

An estimate of this stability criterion can be made by using the 
results of von Kb r man (reference 9). The size of roughness height h 
required just to start to change the skin friction was found to be given 
approximately by 


also 


M- u l 

T-P T ~ 


The agreement between this value and the value obtained by the present 
analysis is good. 

An idea of the accuracy of the boundary- layer model assumed may 
be had by referring to figure 2, which is a replot of the data pre- 
sented in figure 5 of reference 10. The velocity profile shown is that 
in a fully developed turbulent channel flow where the Reynolds number 
based on the half width of the channel is 12,200. The sharp transition 
from laminar profile to turbulent profile is quite evident. If the 

value for — found for the condition of n = 7 is substituted 
into equation (6), there results 
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( 12 ) 


For a value of Re = 12,200, equation (12) yields a value of — - = 0.022, 
o 5 

which agrees well with the experimental value of 0.025 determined from 
the intersection of the straight lines in figure 2. 
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EXTENSION TO COMPRESSIBLE FLOWS 


The extension of the foregoing derivation to compressible flows 
is straightforward. Equation (4) for the ratio of the total stress to 
the laminar stress still holds; however, since the value of v is not 
constant, it must take on its local value dependent on y. The ratio 
of stresses at the edge of the laminar layer is therefore 


so that 


and 



n+1 

! + u o k2& L n 
1 

nVj8 n 



E 


(r - 1) 
k 2 


V L 

UqS 


n 

n+1 


( 13 ) 


^L 

^o 


1 

n(r - 1) V L n+1 
_ k 2 1^8 


(l4) 


When the stress at the wall is evaluated it must be kept in mind 
that the value of the viscosity is not constant through the laminar 
sublayer. The velocity-profile shape assumed is shown in figure 3. 

The stress in the laminar sublayer just below the intersection is again 


assumed to be given by 



which will be taken as 


PL u l/ 6 L* Although 


the viscosity may vary as the wall is approached, the stress must remain 
constant as before; therefore the resulting wall stress is 


1-n 

- = iil^L = ^L u o ( n( r - l) VI 5 * 1 
5 L 5 k 2 UqS 


The skin-friction coefficient becomes 


1-n 2 


n(r - 1) n+1/ v L \n+l p L 

_ k 2 J Wy 


(15) 


c-p- = 2 


(16) 
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When equation (l6) is written in the form 

1-n 2 2 



it may he compared with the value for incompressible flow given by 



If the value of k is not changed appreciably in going to fairly 
high Mach numbers, for example, 0 ^ M ^ or 5, then for the same value 
of Rg and n in each case, 


c f 

Cf M=0 


2 


PL/^y +1 

Wo/ 


( 18 ) 


The magnitude of the reduction in skin friction with Mach number 
as given by equation (17) can be found in the following manner. Since 
the pressure is constant throughout the boundary layer 


PL = ^o 

Po 

and when it is assumed that 

pl /Vr 

Po " Vo) 


equation ( 18 ) becomes 


(19) 


( 20 ) 


n-1 


Cf 


' f M=0 


■ (r ■© 


2_ 

n+1 



n-2m-l 



( 21 ) 
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Since the discussion is for zero heat transfer and since the assump- 
tion of constant total energy in the turbulent layer is not liable to 
introduce too much error, the temperature ratio in equation ( 21 ) may be 
written 



Thus, 


c f 


c 


f M=0 


L 



_ l+ 2 m-n 



( 22 ) 


(23) 


For air where 7 = 1.4, m 


O. 76 , and 


n = 7 > equation ( 23 ) becomes 


c f 


c f 


M=0 



-0.56 


(24) 


Equation (24) shows that the ratio of- skin-friction coefficient for 
the compressible case to the incompressible skin-friction coefficient 
depends principally on Mach number and only to a very limited extent on 
Reynolds number since ul/Uq is not~particularly sensitive to Reynolds 
number. 


The 

equation 


value of_ 



( 13 ) as follows: 


to be put into equation (24) is found from 
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then 



n(r - 1) 
k 2 


2 2(m+l) 

Vo"1 n+1 /TL\ n+1 

^o^] V T o ) 


or 


Thus, 




2(m+l) 

n+1 


n+1 



2(m+l) 


+ AM 2 A^lN 2 
5 Vv 



(25) 


where 



Equation (25) may he solved for (— \ for any Mach number and 

vW 

Reynolds number Rg. These results may then be inserted into equa- 
tion (23) for Cf/c-p 

x i 1 M=0 


Figure 4 is a plot of 



and 


c f / c %=0 


against Mach number 


obtained from equations (24) and (25) for several typical values of Rg 


for a one-seventh-power velocity profile. The insensitivity of skin- 
friction-coefficient ratio to Reynolds number is quite evident. 


It is always interesting for the sake of argument "to see the effect 
on the skin-friction ratio as the Mach number approaches infinity, 
although in this case such an extension is not justifiable. As the Mach 
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number approaches infinity, the term 



approaches unity (see 


equation (25)), but the expression inside the brace in equation ( 24 ) 
approaches l/A, so that 


n-l-2m 


lim 
M— *co 


1 

O 

n(r - l) v q 

\ Cf M= 0 / 

k (l) 2 * * Uo 5 


(n+l)(m+l) 


COMPARISON WITH EXPERIMENT AND OTHER METHODS 


In figure 5 the results of this analysis are compared with the 
experimental results presented by Wilson (reference 7 ) and with the 
experimental and theoretical results of Rubesin, Maydew, and Varga 
(reference 6) as well as with the original suggestion by von Kerman 
(reference ll). The result obtained theoretically by Van Driest (refer- 
ence 8) for a Reynolds number Rjj equal to 7 x 10 & is also compared. 
This result may be compared since, the values of Rg used in figure 5 

cover a range sufficient to encompass the value of 6 that might exist 
at x in the analysis of Van Driest. The agreement between this method 
and the experimental results to date which are described very well by 
the curves of Rubesin and Wilson is good. The curve obtained by the 
method suggested by Monaghan (reference 4 ) almost coincides with the 
experimental curve of Wilson and therefore is in good agreement with 
that of the present method. 


CONCLUSIONS 


A derivation of the form of- the incompressible turbulent skin- 
friction law for an insulated flat plate is made in such a way that it 
may be extended to compressible flows. The following conclusions may 
be drawn: 

(l) The incompressible argument demonstrates that a Reynolds number 

characteristic of the laminar sublayer ujPj^jp. is a constant approxi- 

mately equal to 158. This value may be close to the critical Reynolds 

number of the straight-line velocity profile in the laminar sublayer 

below which all disturbances are damped and above which disturbances 
may feed on the energy of the steady motion. 
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(2) An expression for the ratio of compressible skin-friction 
coefficient to incompressible skin- friction coefficient at the same 
value of Reynolds number based on the boundary- layer thickness on an 
insulated flat plate is derived. The agreement between the compres- 
sible and incompressible skin-friction coefficients found experimentally 
and those predicted by the present method is good. 

(3) The ratio of compressible to incompressible skin friction is 
shown to be rather insensitive to Reynolds number based on the boundary- 
layer thickness and to be principally a function of Mach number. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., February 14, 1952 
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Figure 2 . - Agreement 'between assumed, and experimental profiles from 

reference 10. Rg = 12,200. 








Figure 4.- Plot of 
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from equations (24) and (25). 
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